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PKEFACE 


This  tepcfi  vtras  pcquned  fw  Dr.  Yurii  A.  I^valeako.  Insmute  of  Thcnsopin'sks.  Siberian  Branch 
of  ihe  USS JL  Acaden^’  of  Sdences,  and  Stephen  N.  Flandeis,  Reseatch  Civil  Engineer,  Gvil  and 
Gcoccchnkal  EiKineefuis  Researdi  Branch,  Experimental  Bisinecrins  EMvision,  U^.  Army  Cold 
Regions  Reseanii  and  Engineering  LabcKaioty. 

This  lepod  wiis  originallv*  submitted  to  be  part  of  the  Proceeding  cf the  Workshop  on  In-sttu  Heat 
Flux  Measurements  in  Buildings  that  tvas  held  at  CRREL  on  22  and  23  May  1990.  Its  late  arrival  for 
the  workshop  and  other  consideraiimis  made  publication  as  a  CRREL  Special  Report  appropriate.  The 
report  is  substantially  the  wwk  of  Dr.  Yurii  Kovalenko.  Stephen  Flanders  letvoiked  the  text  in 
cotxretation  with  Dr.  Kovalenko  to  make  it  more  accessibleto  its  readers  inEnglish.  Dr.  Yin-Chao  Yen 
and  Dr.  Virgil  Lunardini  of  CRREL  and  Dr.  Omar  Paiouki  of  the  Queen’s  University  of  Belfast 
provided  technical  reviews. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  purposes.  Gtaiion  of 
brand  names  does  rtot  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial 
products. 
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ILLUSTRATIONS 

Figure 

1.  Determination  of  Cj,  the  ratio  of  contact  surface  area  to  the  total  area  of  the  particle 


suifoces .  7 

2.  Determination  of  a,  the  ratio  of  projected  areas  of  contact  to  the  section  area  of  the 

specimen .  8 


NOMENCLATURE 

a  TbeIIlald^fiils^W^(IIli^^) 

c  ^iecificlieatca|iaci^(J/(lcgK)] 

C  Rdadvebaraze 

C]  ContigiiityCorcoalactdegiee) 

£  Modulus  of  dasticity 

Kyi  Ratio  c^averagedlenipcfaniiegradieots  in 

txnaiy  heterogeneous  system, 
Ki2=lA<r>,|ylA<r>J 

£  Scale  the  proUem(m) 

£e  Lems  number 

Ipr  Process  lengdi  (m) 

Specimen  length  (m) 

/  Molecular  scale  (m) 

ffij  Mass  fraction  ofthe/th  component 

n  Number  of  dimensions 

Nf.  Number  of  particle  contacts 

Np  Number  of  particles 

p  Pressure  of  compaction  (Pa) 

Pcf  Modulus  of  elasticity  or  pressure  of  compac¬ 

tion  (Pa) 

Pf.  Contact  value  .of  modulus  of  elasticiQr  per 

unit  of  contact  area  (Pa) 

9  Local  thermal  flux  density  (W/m^) 

r  Coordinates  (m) 

Sf.  Contact  surface  area  of  a  particle  (m^) 

Scorn  Contact  surface  area  (m^) 

Scompi  Sum  ofparticle  surface  areas  (m^) 

^nom  Nominal  section  area  of  the  specimen  (m^) 

Sp  Surface  area  of  a  particle  (m^) 

Spr  Surface  contact  area  projected  on  the  section 

plane  (m^) 
r  Time  (s) 


T  Tenycrarufe(K) 

BoundaiyvalueQftemperattiieCK) 

T  Heathy  rate  (K/s) 

Vi  Fatial  specific  volume  of  the  ithconqiooent 

(v^/kg) 

z  Distance  normal  to  surfiK%(m) 

a  Rdativecootactsectionofniolded(orsm- 

teroQ  iiiatefiaIs.-The  ratio  of  contact  area 
projected  onto  a  cross  sectioa  to  the  area  of 
the  cross  section  of  die  specimen 
p  Enqmical  parameter,  relating  to  hardness,  p 

» 1.2  forsoftmaterials,  like  aluminum,  and 
P  ~  4  for  hard  metals 

A  Scale  of  a  subre^on  of  the  problem  (m) 

5  Inhomogeneity  scale  (m) 

e  Ratio  ofthe  scale  ofdieinhomc^eneity  to 

the  scale  of  the  problem 

Oj  Angle  between  the /di  contact  plane  and  the 

normal  to  the  section  plane 
V  Ratio  of  thermal  conductivities  for  a  binary 

heterogeneous  system 
n  Porosity 

II<,  Bulk  porosity 

p  Density  (kg/m^) 

Area  or  volume  fraction  ofthe/th  compo¬ 
nent 

X  For  random  In\  and  their  mean  value, 

<lnA>,  X  is  InXr<lnX>,  where  probability  (x) 
is  a  normal  distribution 

\  Contact  thermal  conductivity  (W/m  K) 

A,f  Effective  thermal  conductivity  (W/m  K) 

Maxwell’s  solution  for  effective  thermal 
conductivity  (eq  24) 

<  >  Average 
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llMnnal  Coodiictmly  of  PoixMis  Media  and  Soils 
A'Review  of  Soviet  InTcstigations 

YURn  A.  KOVALENKO  AND  STEFHENN.FLANESRS 


INTRODUCnON 

The  tbennal  conductivi^  of  s(^  and  other  dis¬ 
posed  porous  materials  of  natural  and  aitihcia!  ori^ 
depends,  first  of  all,  on  die  conqxisition  and  geometric 
miciostnictuieofasystem.  Besides  stnicture  and  com- 
porition,  dieimal  conductivity  is  influenced  by  inter- 
phase  (intetcomponent)  interactions.  Heat  transfer  oc¬ 
curs  primarily  by  conduction  in  systems  widi  pores  of 
conventional  size  (<  10~*m),  with  temperatures  rang¬ 
ing  fiom  -SO  to  100°C  and  subject  to  ordinaiy  tem¬ 
perature  gradients.  Gxivective  and  radiative  beat  ex¬ 
change  in  pores  can  be  estimated  (Chudnovskii  1962, 
Dulnev  and  Zarichnyak  1974)  and  is  generally  negligi' 
ble,excq)t  inmaterials  withlaigepores,  such  as  thermal 
insulation,  and  some  stnictural  construction  materials, 
notably  the  large-celled  concrete  found  in  the  Soviet 
Union. 

Porous  media  and  soils  belong  to  a  large  class  of 
Heterogeneous  Media  (HM),  for  which  the  theoiy  of 
heat  transfer  allows  the  investigator  to  take  an  approach 
typical  of  continuous  media:  first,  introduction  of  aver¬ 
aged  characteristics  of  the  medium  and  the  process 
(vdumecontentofcomponents,  temperature, heat  flow, 
heat  capacitance  and  thermal  conductivity),  and  sec¬ 
ond,  formulation  of  principal  physical  laws  in  terms  of 
these  characteristics.  Such  a  continuous  approach,  ap¬ 
plied  independently  of  the  order  or  disorder  of  the  real 
heterogeneous  system,  rests  on  the  assumption  that  the 
scale  of  typical  inhomogeneities  (5)  is  miniscule  com¬ 
pared  with  the  scale  of  the  system  (£,): 

5«L.  (1) 

This  continuous  approach  uses  the  concepts  of  an  infin¬ 
itesimal  volume  (Lifshitz  and  Pitaevskii  1979)  and  of 
incorporated  continua.  Any  such  averaged  description, 
as  an  alternative  to  a  detailed  description  of  processes  in 
HM,  considerably  simplifies  the  problem  solution,  pro¬ 
vided  that  the  regular  microstmcture  is  known.  It  allows 
the  solution  of  practical  problems  when  the  inner  struc¬ 


ture  is  non-regular  and  unknown  (Ity  the  statistic^ 
mediod  in  Landau  and  Li&hitz  [1976]).  The  difficulty 
of  describing  processes  in  this  case  is  in  obtaining 
averaged  equation  coefficients  for  energy,  momentum 
and  mass  (analogous  to  obtaining  closure  for  hydrody¬ 
namic  equations  [lifshitzand  Pitaevskii  1979]).  Bang 
random  variables  fn-  5  ^L,  these  coefficients  become 
single-valued  characteristics  of  the  medium  when  the 
conditions  of  eq  1  are  satisfied  (the  so-called  self- 
aveia^g  of  specific  extensive  quantities  and  kinetic 
coefficients  as  5/L  0  fiom  Lifshitz  et  al.  [1982]). 

Several  fiuther  concepts  of  HM  classificatimi  are 
useful.  Many  classifications  by  different  indicators  ex¬ 
ist  (for  example.  Chudnovskii  1962,  Dulnev  and 
Zarichnyak  1974,  Dulnev  1979,  Heifitz  and  Neimark 
1982).  The  principal  indicators  are: 

1.  The  occurrence  of  binary  or  multi-component 
mixtures. 

2.  The  occurrence  of  inert  mechanical  mixtures, 
alloys,  solutions,  HM  with  chemical  interactions,  etc. 

3.  The  topology  of  the  microstructure — including 
the  degree  of  surface  contact  and  the  distribution  of 
component  sizes  within  the  medium,  ranging  frcmi 
uniform  to  varied. 

4.  The  degree  of  order  or  randomness  in  the  HM. 

5.  The  uniformity  or  variety  of  the  components  with¬ 
in  the  HM,  ranging  from  weakly  to  strongly  heteroge¬ 
neous.  This  includes  the  physical  qualities  of  the  com¬ 
ponents  and  their  arrangement  within  the  medium. 

These  classifications  are  conventional  and  in  prac¬ 
tice  the  given  types  of  systems  are  mbced,  with  a  high 
degree  of  disorder  of  components  being  typical.  In  this 
respeetporous  media  andsoilsaresimilartoothertypes 
ofHM. 

The  first  investigations  in  which  the  continuous  ap¬ 
proach  was  extended  from  homogeneous  gases,  solu¬ 
tions  and  liquids  to  heterogeneous  media  with  inner 
structure  are  in  works  Ity  Maxwell  ( 1 873)  and  Rayleigh 
(1892),  which  focused  on  calculation  of  the  effective 
electric  field  in  a  medium  with  spherical  inclusions  that 
are  located  at  the  nodes  of  a  cubic  lattice.  Burger* 


(4) 


studieddiecaseofellipsoidalincIuaons.Ihl887Ar- 
liienius*  studied  the  aoalogotisproUein  for  the  vucbs- 
iiyofUnaiy  liquid  solutioos  and  suggested  tbefonnula 
for  the  "logarithmic  law  ofviscosilyniixing.’’Later  the 
same  approach  appeared  again  in  Uchtenecker’s  woik 
(1909,  1924,  1926,  1929)  on  conductivity.  Eydcea* 
(1912)  q)paiently  was  the  first  to  study  the  equivalent 
protdem  for  heat  transfer.  Analogous  proUans  coo- 
nected  with  various  physical  fields  were  alsostudiedby 
Wiener*  (1912),  Fricke*  (1924),  Voight*  (1928),  Re- 
uss*(1929),Bmggeman(193S)  andothers(Oddevskii 
i951,  Chudnovskii  1962,  Vasiliev  andTanaeva  1971, 
DulnevandZarichnyak  1974,DulnevetaL  1976,Bakh- 
valoy  and  Panasenko  1984,  Shvidler  1985).  A  modem 
introduction  to  the  formal  study  of  physical  character¬ 
istics  that  are  analogous  to  conductivity  (the  dielectric 
constant  and  magnetic  permeability,  viscosi^,  electri¬ 
cal  and  thermal  conductivity,  difhisicm  and  filtraticm 
coefricients,  and  moduli  of  elasticity),  is  found  in 
Odelevskii  (1951). 

Investigations  conducted  in  the  U.S.S.R.  may  be  di¬ 
vided  into  two  large  groups:  1)  ad-hoc  experiments  in¬ 
vestigating  the  thermal  conductivity  of  specific  materi¬ 
als  andmedia,and2)theoreticalcalcuIationstodeteimine 
general  solutions  for  effect!  ve  thennal  conductivity  and 
overall  conductance  in  wide  classes  of  HM. 


CONTINUOUS  DESCRIPTION  OFTHERMAL 
PROCESSES  IN  HETEROGENEOUS  MEDIA 
AND  THE  CONCEPT  OF  EFFECTIVE 
CHARACTERISTICS 

The  fundamental  distinction  between  averaged  de¬ 
scriptions  of  HM,  and  macroscopic  descriptions  of  mo¬ 
lecular  systems  rests  with  the  heterogeneity  of  the  HM 
being  much  larger  than  molecular  scale  (/): 

5»/.  (2) 

The  averaging  technique  used  allows  one  to  approxi¬ 
mate  a  continuous  medium  with  the  same  form  of  con¬ 
servation  of  energy  equations  as  for  the  HM  and  its 
internal  heterogeneities,  each  having  its  own  local 
physical  coefficients. 

Even  if  the  contact  area  on  the  interface  between 
components  in  a  micro-heterogeneous  system  is  large, 
one  may  often  neglect  the  interaction  of  components 
within  it  over  a  large  temperature  span.  In  such  temper¬ 
ature  intervals,  HM  are  mechanical  mixtures  and  heat 
transfer  is  described  by  the  pure  conduction  equations: 

div^-Hcp  — =  0  (3) 

_  dt 

*  Citation  not  available. 


q=-XffddT 


(5) 

Hero  qxcific  heat  c(r),  density  p(r)  and  tbennal 
conductivity  have  n^dly  changing  coordinate 
fiii^ons  (with  periodical  or  random  variables).  As 
cp(^  and^^arodiscontinuousinHM.tbesolutionof 
the  eq  3  tbrou^  5  is  considered  to  be  a  general  solution 
satisfyii^  corresponding  integral  relations  (Bakhvalov 
and  Panasenko  1984).  It  is  necessary  to  assume  a 
temperature  and  heat  flux  continuum  along  interfaces 
between  components 

[T]=0. 


where  []  isthefiinctionjumpon  thesuifaceofacontact 
and  dr/dz  is  the  partial  derivative  of  temperature  nor¬ 
mal  to  the  surface.. 

hi  priiKiple,  a  detailed  knowledge  of  the  HM  micro- 
stmcturo  that  deteiminesthe  functions  c(r),p(r)  andA(r) 
would  allow  one  to  solve  the  problem  accurately.  With 
the  assumption  of  eq  1  and  if  the  structure  is  rather 
complicated,  then  such  a  detailed  solution  would  be 
difficult  to  model  on  a  computer.  However,  such  a 
computation  would  rely  on  a  detailed  knowledge  of  the 
microstructure  of  the  HM,  which  is  unknown  in  detail. 
Instead,  collective  general  properties  of  the  microstnic- 
ture  are  typically  known.  Thanks  to  this  knowledge,  a 
detailed  solution  isnotusuallyrequiredandanaveraged 
description  of  the  medium  and  the  processes  taking 
place  within  is  typically  the  only  appropriate  method.  It 
is  also  advisable  to  average  the  processes  in  simple, 
periodic  media,  as  was  done  in  Bakhvalov  and  Panasen¬ 
ko  (1984)  and  Dulnev  et  al.  (1976). 

Now  we’ll  cover  isotropic  and  macroscopically  ho¬ 
mogeneous  systems,  that  is,  the  systems  in  which  the 
scale  of  heterogeneity  5  and  the  averaged  parameters 
(thevolumcofthecomponentsandsoon)donotdepend 
on  global  coordinates.  In  terms  of  the  multiple  scales 
method  that  is  applied  to  averaged  processes  in  period¬ 
ical  media  (Bakhvalov  and  Panasenko  1984),  these  so- 
called“slowvariables”areonthescaleofL,unlike“quick 
variables,"  which  are  inside  heterogeneities  and  are  on 
the  scale  of  8. 

Averaging  can  be  carried  out  by  several  methods. 
The  method  strongly  depends  on  the  microstructure  of 
the  HM.  It  may  be  1)  statistical  averaging  over  an 
ensemble  of  similar  systems,  where  cp(r)  and  X{r)  are 
random  functions,  2)  averaging  over  a  physically  small 
volume  A^,  8  «  A  «  L,  or  3)  averaging  over  an  ele¬ 
mentary  cell  that  is  in  a  periodic  structure.  No  matter 
how  the  averaging  is  done,  it  essentially  uses  the  disin¬ 
tegration  of  the  solution  from  Bakhvalov  and  Panasen- 
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ko  (1984)  or  tile  pz.aboIic  opeiator  of  equations  from 
Koz]ov(1978).ShvidIer(1986)andZhikovetaI.(1981). 
sequentially,  into  a  power  series  of  the  parameter  e  = 
c/L.  Averaged  equations  of  the  zeroth  order  corre¬ 
sponding  to  eq  3  and  4  are  expressed  as  follows 


div<^  <cp>  §53^=0 
dt 

(6) 

<q> = -Xgf  grad  <r> 

(7) 

where  the  operations  div  and  grad  are  taken  only  by 
global  variables.  In  eq  6-8  the  symbols  <  >  mean 
averaging  1)  over  an  ensemble  of  similar  systems  with 
randcnn  fields,  cp(r)  and  A(r),  2)  of  a  physically  small 
volume  A^  or  3)  using  the  periodicities  of  regulaiities 
within  the  structure.  In  order  to  solve  for  <7^  approxi¬ 
mately  in  eq  6  and  7  to  solve  for  T  in  eq  3-5 

max  I T -  <r>l  <  CqE 

mustbe  satisfied,  whereCoisaquantityofthefirst  order 
of  magnitude.  One  must  establish  boundary  conditions 
so  that  they  coincide  with  their  averaged  analogue 

<r>lB  =  r(rB./).  (8) 


We  now  look  at  the  concept  of  effective  thermal 
conductivity  X«f,  which  is  discussed  as  a  special-case 
ancillary  problem  in  Bakhvalov  and  Panasenko  (1984), 
Shvidler  (1985)  and  Shvidler  (1986),  Characterizing 
is  not  trivial  and  represents  the  main  challenge  in  this 
whole  topic.  Many  works  are  devoted  to  the  theoretical 
determination  of  However,  experimentation  is  the 
most  reliable  way  for  determining  X^f. 

The  other  thermal  coefficient  derived  from  averaged 
equations,  average  theimal  capacitance  <cp>  in  eq  6,  is 
determined  simply  by  additivity.  In  fact,  by  dividing 
regions  of  integration  into  separate  subregions  occu¬ 
pied  by  individual  components,  it  is  easy  to  show  that 


where  c\  and  pci  =  specific  heat  capacity  and  density 
/  =  component 

(!)i  =  volume  fraction /-component  (aver¬ 
aged  over  the  characteristic  volume 

a3). 


By  introducing  partial  density  (the  concentration  of  the 
/‘•’-component),  pi = (jiiPci,  equatim  9  may  be  expressed 
as 


<cp>  =  X‘^iPi- 

Anotherderivation  of  averaged  volume  heat  capacity  is 
possible  by  using  the  effective  specific  heat  and  an 
averaged  density  of  the  HM.  Let’s  first  determine  the 
mass  content  of  the  components  in  their  mixture 

i 

and  then  the  partial  specific  volume 

Vi  =  mi/pci. 

By  using  vj,  we  determine  the  volume  content  of  the  i’'’- 
component  in  a  manner  analogous  to  eq  1 1 

<l>i  =  Vi/X'’i- 
i 

Then,  introducing  the  effective  specific  thermal  capac¬ 
ity  of  the  HM  as 

Ccf  =  Xw  (12) 

i 

and  the  average  density 

<p>=pfe=SPi=(2'’iJ-'  (13) 

from  eq  10-13  it  is  easy  to  show  that 

<cp>  =  Cef<p>-  (14) 

Note  that  the  effective  specific  heat  Qf  does  not  coin¬ 
cide  with  the  averaged  specific  heat,  as  determined  by 
the  mean  of  individual  specific  heats,  namely  Ccf  <c>. 
The  same  is  tme  for  Xgf  <A>.  In  the  particular  case  of 
layered  HM,  as  we’ll  see  below,  =  <X>.  These 
quantities  agree  only  for  volume-specific  properties. 

Note  the  following  limitation  on  the  applicability  of 
eq  6  and  7  to  non-steady-state  processes.  Their  validity 
may  be  accepted  only  when  the  local  temperature  within 
the  medium  differs  only  a  little  from  the  temperature  of 
the  local  thermal  equilibrium.  They  are  unconditionally 
valid  only  when  the  medium  is  close  to  steady  state,  that 
is,  at  dT/dt = 0  in  Rubinshtein  ( 1 948).  Hence,  it  follows 
that  they  can  describe  only  rather  slow  non-steady-state 
thermal  processes  with  typical  response  times  much 
longer  than  the  time  constants  of  the  heterogeneities 
within  the  medium,  that  is 

(15) 

^in 

where  Amin = min  (a,)  and  can  be  described  by  averaged 
eq  6  and  7.  The  condition  (eq  15)  can  also  be  obtained 
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from  the  estimate  of  the  teims  of  eq  6  and  7  using  the 
inequality  of  eq  I  (Buevich  1973).  The  effective  ther¬ 
mal  conductivity  in  eq7  for  non-steady-state  processes, 
tested  against  the  condition  in  eq  1 S,  coincides  with  the 
steady-state  effective  thennal  conductivity. 

If  one  cannot  adhere  to  the  assumptions  of  local 
'  thennal  equilibrium,  then  one  has  to  resort  to  a  more 
detailed  description — forexampIe,amulti-temperature 
model  (Rubinshtein  1948,  Schvidler  1986).  In  such  a 
model,  heat  transfer  terms  that  are  linear  functions  of 
temperature  are  added  to  eq  6  for  each  temperature 
within  the  individual  components. 

Summing  up,  let  us  reformulate  the  conditions  of  the 
applicability  of  homogenization  of  both  steady-state 
and  non-steady-state  thermal  problems  in  HM,  including 
mechanical  mixtures  and  also  dispersed  systems  with 
possiblephaseandchemical  reactions  (Kovalenko  1 987). 

Withtheinequalityofeq  1  as  the  principal  condition, 
the  smallest  of  the  following  two  dimensions  must  be 
chosen  as  the  typical  problem  scale:  1 )  minimal  size  of 
specimen  Z-jp,  or  2)  the  typical  process  length  Lpr,  that 
is:  L  =  min  (Ljp,  Lpr).  TTie  condition  L»8  must  be 
observedforanysteady-state(thisusuallyallowsL=Lsp) 
or  non-steady-state  process.  In  this  latter  case  the  con¬ 
dition  of  eq  1  can  be  rewritten  by  the  substitution 

^  ^  ^char'^ 

as  the  condition  of  eq  15,  in  which  the  L  =  Ijp.  the 
characteristic  time  is  simply  equal  to  Lf^/a.  When  L  = 
Lpr,  the  characteristic  time  can  be  expressed  by  other 
parameters.  So,  forHM  undergoingphase  endothermal 
change  from  thawing,  =  AT/T  can  be  introduced, 
where  Aristheliquidcharacteristictemperaturediffer- 
ence  in  the  phase  diagram  (melting  phase),  and  where  T 
is  the  heating  rate.  The  value  /char  "Ot  only  represents  eq 
15  but  exceeds  the  diffusion  time  in  the  heterogeneity 
scale.  The  latter  condition  is  thanks  to  the  fact  that  a 
Lewis  number,  Le  <  10“^  «  1,  in  condensed  media 
appears  sufficient  for  homogenized  modeling  of  the 
medium  (in  this  case  the  condition  of  eq  15  is  deliber¬ 
ately  observed).  With  such  a  homogenization,  we  un¬ 
derstand  the  description  of  the  medium  with  the  help  of 
lumped,  effective  characteristics. 

THEORETICAL  METHODS  FOR  THE 
DETERMINATION  OF  EFFECTIVE 
THERMAL  CONDUCTIVITY 

By  definition  the  effective  thermal  conductivity  is  a 
coefficient  that  relates  the  averaged  values  for  heat  flux 
and  the  gradient  of  the  average  temperature  (eq  7) 

<^>  =  -Xefgrad  <T>. 


Most  calculation  methods  are  based  on  this  defini¬ 
tion.  The  problem  then  is  the  determination  of  the  fields 
of<^>i  and  <T>i  averaged  by  the  components  within  the 
HM. 

Exact  solutions  and  approximate  methods 

Exact  solutions  to  the  problem  of  the  definition  of  the 
effective  theimal  conductivity  are  known  only  for  a  few 
cases.  In  the  first  place,  these  are  .one-dimensional 
layered  structures  for  which  the  longitudinal  conduc¬ 
tivity  is  as  follows 

X||=<X>  (16) 

and  transverse 

(17) 

For  more  complex  structures,  especially  for  non¬ 
regular  ones,  it  is  complicated  enough  to  calculate  the 
effective  thermal  conductivity  on  the  basis  of  the  infor¬ 
mation  about  field  structures  of  <q>i  and  <T>i,  so  exact 
solutions  are  very  rare.  However,  one  may  occasionally 
avoid  awkward  and  sometimes  unsuccessful  calcula¬ 
tions  for  systems  that  are  not  strongly  heterogeneous. 
This  possibility  relies  on  using  variational  principles  for 
estimating  the  bounds  of  effective  characteristics 
(Dykhne  1967,  Shvidler  1985).  The  simplest  expres¬ 
sion  of  such  bounds  is 

(18) 

which  was  probably  determined  for  the  first  time  by 
Wiener  and  later  was  ascertained  by  Hill*  ( 1 964)  for  the 
modulus  of  elasticity  and  pliability  and  by  Dykhne 
(1967)  for  conductivity. 

Dykhne  (1970)  obtained  an  interesting  result  in  his 
work,  in  which  he  proved  that  fora  flat  isotropic  system, 
covered  on  the  average  by  geometrically  equivalent 
regions  with  different  conductivities  A,]  and  A,2,  the 
effective  conductivity  satisfies  the  functional  equation 

^f(<l>i)  •  (<t>2)  =  ^1^2 

where  (j);  is  the  area  fraction  (in  a  two-dimensional 
model)  for  regions  with  conductivity  Xj  and  <}>i+<l)2  =  1. 
Hence  at  (!)i  =  (i)2  =  V2  follows 

(19) 

In  the  above  mentioned  work  by  Dykhne  and  also  a 
paper  by  Kozlov  (1979),  the  case  when,  for  the  random 
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variable  X,  the  distribution  of  the  valuex = lnA,-<lnA> 
is  an  even  function  ofx,  is  studied  fordi^eient  possibil¬ 
ities.  The  precise  value  of  tlie  effective  thermal  conduc¬ 
tivity  is  obtained  as 

Aef  =  exp  <lnX>.  (20) 

This  result  with  its  particular  case  of  eq  19  represents  the 
empirical  logarithm-mixing  formulaby  Arrhenius  (also 
suggested  later  by  Lichtenecker)  which,  thus,  is  accu¬ 
rate  for  two-dimensional  bicomponent  systems  with 
specialized  conductivity  distribution  functions.  Equa¬ 
tion  20  can  be  represented  by  (Shvidler  1985) 

Aef=(<X><A,->>->)'/2.  (21) 

This  is  a  geometric  mean  that  is  between  the  arithmet’c 
and  harmonic  means  for  thermal  conductivity. 

In  the  work  in  Ivanov  ( 1 967),  eq  21  is  generalized  to 
other  dimensions  /t  with  the  help  of  the  additional 
results  of  the  perturbation  method 

.  (22) 

At  n  =  1  this  fomiula  is  accurate  for  any  distribution,  at 
«  =  2  it  is  accurate  if  the  distribution  of  %  is  even  and  at 
« = 3,  according  to  Shvidler  (1985),  “one  can  expect  the 
formula  to  be  accurate  enough  for  normal  distributions 
oflnX.” 

The  approach  made  in  the  early  works  of  Maxwell 
(1873)  and  Rayleigh  (1892),  which  present  an  approx¬ 
imation  of  a  small  concentration  of  inclusions  dispersed 
either  regularly  or  randomly  within  the  matrix,  gave  rise 
to  a  great  number  of  publications  in  which  different 
asymptotic  formulae  for  were  obtained.  Maxwell 
(1873)  derived  the  first  term  of  a  decomposition  of 
effective  thermal  conductivity  as  a  function  of  the 
volumetric  fractions  of  the  inclusions 


v(1-<)),)  +  2-k)), 


(23) 


where  the  subscript  1  pjfers  to  the  inclusions,  2  refers  to 
the  matrix,  and  v  -  X\fk2-  Rayleigh  (1892)  studied  the 
problem  of  spherical  inclusions  at  the  nodes  of  a  regular 
cubic  lattice.  He  obtained  the  first  two  terms  of  an 
expansion  of  X^fby  degrees  of  (l»i;  the  first  and  principal 
term  of  this  expansion  disregarded  the  mutual  influence 
of  inclusions,  and  coincided  with  Maxwell’s  (1873) 
solution  (eq  23).  Rayleigh’s  solution  was  rendered 
more  precise  by  adding  two  more  terms  of  the  expan¬ 
sion  in  Berdichevskii  (1979). 

Hashin  and  Strikman  (1962)  narrowed  the  limits  of 
effective  parameters  beyond  those  in  eq  18  using  the 


variational  principle.  For  the  isotropic  two-phase  sys¬ 
tem  with  ^  X2,  the  limits  coincide  with  asymptotic 
formulae  by  Maxwell  (1873) 

^<l>l+2+  4>2 

/•I 

^(H24.i-)+2<t>2 

=  5^2  *^2 - .  (24) 

^<t>2  +  2+<t>i 
4-2 

This  is  not  a  casual  coincidence.  In  fact,  as  shown  in 
Berdichevskii  (1979),  if  the  thermal  conductivity  of 
inclusions  X]  is  less  than  that  of  the  medium  X2,  then  the 
accurate  value  of  effective  conductivity  of  a  periodical¬ 
ly  structured  medium  with  spherical  inclusions  is  less 
than  or  equal  to  the  first  term,  of  the  Rayleigh- 
Berdichevskii  formula  in  Berdichevskii  (1979)  that 
coincides  with  Maxwell’s  solution  (eq  23).  If  Xj  S  X2, 
then  Xgf  ^  Xjjj. 

Unlike  the  universal  bounds  forA^f  in  eq  1 8,  Hashin- 
Strikman’s  expression  (eq  24)  forthe  three-dimension¬ 
al  heterogeneous  system  gives  approximate  bounds, 
which,  however,  are  considerably  narrower  than  the 
former  (eq  18).  Using  such  bounds  it  is  possible,  when 
necessary,  to  constmct  approximate  solutions,  as  in  eq 
24.  Shvidler  (1985)  further  recommends  that  if  compo¬ 
nent  connectivity  is  identical  (mutually  penetrative  com¬ 
ponents),  then  we  can  take  as  an  approximate  value  of 

^cf=(^m+^iTi)/2. 

For  isotropic  matrix  HM,  he  recommends  the  assump- 
tionXjf  =  X^,  if  thematrix  thermal  conductivity  is  more 
than  the  inclusion  conductivity,  and  X^f  =  X“  in  the 
opposite  case.  If  no  single  component  is  simply  con¬ 
nected,  but  the  components’  connectivities  are  differ¬ 
ent,  then  the  estimate  for  X(.f  is  X^f  =  Vx^Xj^  ■ 

The  perturbation  method  for  determining  thermal 
conductivity  in  Landau  and  Lifshitz  (1982)  is  similarto 
the  Maxwell  (1873)  method  for  approximating  small 
concentrations  of  inclusions.  The  perturbation  method 
gives  good  results  for  any  concentration  of  components 
but  only  for  weakly  heterogeneous  systems. 

In  addition  to  the  perturbation  method,  the  self- 
consistent  effective  field  method  is  another  important 
technique.  In  a  number  of  cases  it  gives  better  results. 
The  first  self-consistent  parameters  were  probably  cal- 
culatedby  Bruggeman  ( 1935).  Later,  Odelevskii  (1951) 
in  the  U.S.S.R.,  Landauer  (1952)  and  others  developed 
this  method  for  other  systems.  The  efficiency  of  the  self- 
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consistent  field  method  is  shown  by  direct  numerical 
computations  (Shvidler  1983)  and  comparison  with 
experiments  (Odelevskii  1951).  Shvidler  (1985)  gives 
tables  for  the  values  of  AefAi  calculated  by  this  method 
for  different  values  of  X2A1  and  <|)i  and  for  different 
dimensional  ratios  of  spheroidal  inclusions. 

Note  that  in  the  strongly  heterogeneous  HM,  perco¬ 
lation  causes  threshold  ejects  (Efros  1982,  Shvidler 
1985).  One  can  describe  conductance  behavior  near  the 
threshold  in  such  systems,  using  power  laws  for  scaling 
according  to  percolation  theory.  Critical  indexes  and 
values  of  the  threshold  concentration  for  simple 
problems  can  be  obtained  analytically,  but  more  often 
they  are  computed  by  the  Monte  Carlo  method.  Self- 
consistent  parameters  at  v  =  0  or  at  v = «» in  some  con¬ 
centration  areas  get  non-physical,  negative  values.  While 
this  is  certainly  the  result  of  approximation,  it  is  also  an 
indicator  for  the  percolation  effect  (Shvidler  1985). 
Thus,  the  formula  for  binary  HM  with  equal  compo¬ 
nents  derived  in  Odelevskii  (1951)  by  the  effective 
media  method 

>-ef  =  ->{(3<l>rl)+v(3.1)2-l)  + 

4 

V[(3(|)2-1)+ v(3(t)2-l)p+8v)  (25) 

where  v  =  A,2A.i  gives  the  estimate  of  percolation 
threshold  (jrc  =  1/3  as  v  0 

when  (j)i  S  (t)c  Xef  =  0 

when  <1)1  ><!)c  Xef  =  X,i(3(j)i-l)/2,  (26) 

In  engineering  practice  the  methods  of  structural 
modeling  (Chudnovskii  1962,  Vasiliev  and  Tanaeva 
1 97 1 ,  Dulnev  and  Zarichnyak  1 974,  Dulnev  et  al.  1 976, 
Dulnev  1979)  are  often  used.  The  real  HM  structure  is 
modeled  by  the  most  suitable  ordered  structure  in  which 
an  elementary  cell  is  separated,  and  its  thermal  conduc¬ 
tivity  is  calculated  precisely  using  a  computer  as  in 
Bakhvalov  and  Panasenko  (1984)  or  approximately  as 
in  Dulnev  and  Zarichnyak  (1974).  Approximate  calcu¬ 
lations  of  Xef  for  the  elementary  cell  are  usually  made 
by  dividing  a  cell  into  separate  fragments,  containing 
only  one  of  the  components  with  surfaces  that  are  ortho¬ 
gonal  to  the  direction  of  heat  flow.  Then  the  thermal 
resistance  of  acomplete  elementary  cell  is  made  of  par¬ 
allel  and  serial  connections  of  fragment  resistances  and 
calculated  by  the  rules  of  Kirchhoff’s  chain.  A  number 
of  formulae  sufficiently  useful  for  engineering  calcula¬ 
tions  are  obtained  by  this  method  in  Dulnev  and  Zarich- 
ny  (1974).  Thus,  a  model  was  suggested  for  structures 
with  interpenetrating  components,  in  which  one  of  the 
components  presents  a  cubic  lattice  or  frame  with  bars 


of  a  constant  square  section.  To  determine  the  conduc¬ 
tivity  of  such  a  system  by  the  method  described,  the 
estimates  of  \{  were  obtained  with  the  help  of  two 
different  fragment  thermal  resistance  connections  with¬ 
in  the  elementary  cell 


\^=7ilc\v(l-cf  +  2vC(l-CXvC+l-C)“'] 

r 


^isot  “^1 


1-C 


C^+v(l-C^)  C(2-C}fv(l-C^)j 


where  the  subscripts  ad  and  isot  refer  to  adiabatic  and 
isothermal,  respectively,  and  C  is  the  relative  bar  size, 
which  is  a  root  of  the  equation 


2C3-3C2+1-(|)2=0. 


With  the  arithmetic  mean  of  these  values  from  Dul¬ 
nev  and  Zarichnyak  (1974),  one  can  determine  the 
effective  value  of  k  for  any  relationship  of  the  compo¬ 
nent  thermal  conductivities,  the  result  being  rather  near 
to  the  results  of  numerical  computations.  Maximum 
calculation  error  is  attained  when  v  =  0  or  v  =  «>,  yet  it 
doesn’t  exceed  -15%  (Dulnev  and  Zarichnyak  1974). 

These  methods  include  the  averaged  element  tech¬ 
nique  (Dulnev  and  Zarichnyak  1974).  This  technique 
introduces  a  set  of  assumptions  about  the  HM  structure 
and  then  averages  the  parameters  of  the  structure.  The 
complexity  and  subjectivity  of  the  procedure  to  define 
such  an  averaged  element  substantially  explains  why 
this  method  was  not  widely  adopted. 

In  some  applied  problems,  it  is  necessary  to  take  into 
account  the  statistical  distribution  of  heterogeneities  in 
HM  leading  to  deviations  in  thermal  conductivity  val¬ 
ues  from  the  average  value  A^f.  Zhirov’s  work  (1976) 
focuses  on  this  problem. 

The  success  of  the  application  of  all  these  formulas 
to  real  HM  depends  on  how  well  one  has  characterized 
the  system  structure,  e.g.,  the  distribution  of  A  for  ran¬ 
dom  mixtures  or  for  geometrical  structures.  Variational 
estimates  based  on  universally  used  principles  are  rela¬ 
tively  simple  and  require  minimal  information  about  the 
HM  studied.  Therefore,  they  apply  to  wide  classes  of 
HM.  This  is  their  advantage  and  at  the  same  time  a 
limitation  on  their  efficiency  when  they  lead  to  wide 
bounds  in  estimating  effective  conductivity.  One  can 
obtain  positive  results  provided  that  0.2  <  v  <  5  for  the 
Hashin-Strikman  estimate  and  0.5  <  v  <  2  for  the  uni¬ 
versal  relationship  are  both  true.  Attempts  to  include  in 
the  analysis  more  detailed  information  about  specific 
features  of  the  systems  studied,  in  order  to  narrow  the 
limits  of  estimates,  instead  make  the  variational  esti¬ 
mates  method  considerably  more  cumbersome. 
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|V<r>|  ' 


(27) 


where  V  =  "hifkx,  and  A’i2= !  V<T>iM  V<r>j  is  the 
rado  of  the  averaged  temperature  gradients  in  die  two 
components  of  a  binary  heterogeneous  system  (after 
Shvidler  1986).  A  number  of  prt^ietties  of  ^{2 can  be 
set  Thus,  using  eq  18  it  is  easy  to  show  that  this  value 
is  limitedby  the  boundsforAri20fl,vand,  hence,  when 
V  =  1,  then  A'i2  =  1.  With  die  help  of  such  boundary 
relationships  for  dielectricpolarizadonfincHn  Odelevskii 
(1951)  and  Landau  and  Lifshitz  (1982),  we  obtain  the 
value  of  the  derivative  when  v  =  1 

§^11 =i. 

9v  3 

There  exist  generalized  experimental  data  from  Malter 
et  al.  (1980)  that  indicate  that  for  a  wide  class  of 
systems,  when  V = 0  or  V = oo  the  boundary  values  of Tfi  2 
(v)  do  not  depend  on  the  volume  fraction  over  a  wide 
range  and  satisfy  the  condition 

0.4^/i:,2(0)  +  fi:,2{«)<0.67 

U,2=^21')- 

Assuming  similar  properties  of  the  ratio  /(n,  Malter 
et  al.  (1980)  suggested  a  four-parameter  approximate 
equation  for  fifi2  (v),  considering  it  independent  of  the 
concentration.  Assuming  that  Kyi  is  independent  of  (|ij, 
it  is  possible  to  determine  this  constant  for  the  given  v 
and  from  eq  27  if  the  value  of  has  been  measured  for 
certain  values  of  ({)j;  then  for  the  other  values  of  (|ii  the 
tlrerraal  conductivity  can  be  calculated  by  using  the 
derived  value  ofAT^. 

Alsoof  interest  are  the  works  of  Serykh  and  Kolesni¬ 
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kov  (!9S2)  and  KoSesrukov  (1935).  witkii  SE^gesi  Gse 
of  stereologscal  asxdysEs  fos*  obcahtias  the  aa^ioeg! 
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(28) 

HereC]  iscalkdcontiguityord^ieeofcositacL  Iiistlg 
ratio  of  the  contact  surface  area  to  the  total  area  of  the 
particle  surface  (Fig.  1) 

Cx  =  2Sc^S.^  (28a) 

Nc 

whcrcSeaa= 

^=1 

5c  =  contact  surface  area 
Nc  =  number  of  contacts 
and  A'p 

5p=  p^icle  surface  area 
A^p=  number  of  particles. 

C|  is  determined  from  the  stereological  analysis  of  a 
microseciton  specimen.  This  method  can  sometimes 
simplify  thermal  experiments  for  porous  granular  sys¬ 
tems,  but  usually  direct  measurement  of  thermal  con¬ 
ductivity  is  simpler. 

In  Balshin’s  work  ( 1 948, 1 972)  the  followmg  formu¬ 
la  for  the  conducti  V  ity  properties  of  porous  matenals  is 
suggested 

(29) 

in  which 

~  -pi  '  5^nom  (29a) 


(28b) 


(28c) 
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Pc{=P<!a-  (30) 

HcrcPefniay  betbemodulusofelasiidlQ  £.tbesiiengih. 
or  the  piessuie  of  c(Mnpacti(Hi p.  Each  of  these  repmseni 
a  resistance  prc^ity  for  forming  the  material.  Fc  is  the 
corresponding  contact  value  for  this  property.  For  £,Fc 
IS  the  corresponding  modulusof  elasticity  of  the  frame* 
matenal.  For p,  F^  is  acorresponding  value  between  the 
matenal  frame  strength  and  the  yield  point  of  this 
matenal  (Babhin  1948. 1972). The  weight  rftheexper* 
imental  evidence  in  Balshin  {, 1 972)  suggests  an  approx¬ 
imation  of  the  contact  section  value  through  measured 
parameters  either  for  porosity  11  and  the  initial  bulk 
porosity  rio  or  corresponding  values  of  the  volume 
fraction  of  the  solid  phase 

a=(i-n)2(i-n/no)P  od 

where  p  depends  empirically  on  the  hardnesses  of  the 
material  species. 

The  stereological  definition  of  the  contact  section  is 


*  The  fnnie  refers  lo  ihe  solid  skeleton  of  a  porous  material,  exclusiv  e 
of  the  voids  within  it. 


(iesmbed  fov  eq  2^  Us^  this  deiciidoa.  Sesy-kh  atad 
Kolesskcr  (1932)  aad  KoSeatsEkov  (1985)  oicaected 
the  values  Cj  sad  or  the  of  dae  stxceoSo^cal 

agpKgidi 

C,  =  c?fd,.  (32) 

Hezaoe.  ax  can  6smc  the  fonntrSa  for  the  effective 
mennal  coenimhiiy'  idetakal  to  eq  29  from  eq  2S. 
This.  OEK  riti^  march  both  approaches  and  desimiraf 
aby  steneolo^cal  mears. 

The  formulas  fcK^  7.^1  in  eq  28  and  29  are  based  on 
some  lax  assumptiocis  that  tSo  txM  aoccuni  fw  the  phys> 
teal  state  of  the  contacts.  This  pmbably  overestimates 
ihetheimalaKtdnctivity’insomecases.fOTexample.for 
aKn{Kessed  metal  powders  with  high  pmoriiy.  These 
systems  were  studied  by  Aleksarxlrov  ct  al  (1985). 
Gruzdev  and  Kovalenko  (1987, 198$a)  and  Gruzdev  el 
aL  (1989),  using  crxnpressed  pow  ders  of  nickel,  zirco¬ 
nium  and  their  mixtures  with  aluminum  powder  as 
examples.  In  .Aleksandrov  et  al.  (1985)  and  Gruzdev  el 
al.  (1989).  experiments  determined  that  a  good  analogy 
exists  between  heat  transfer  and  mechanical  stress 
transfer  in  these  systems,  as  deierm  ined  by  the  degree  of 
contact  present.  A  lattice  model  of  a  powder  body 
demonstrates  the  analogy  theoretically  in  Gruzdev  and 
Kovalenko  (1987,  1988a). The  thermal  versus  mechan 
ical  analogy  betw  een  the  molding  pressure  and  thermal 
conductivity  is  based  on  the  similarity  of  dependence  on 
specimen  porosity.  Employing  the  relationship  of  po¬ 
rosity  tomcchanical  properties  in  cq  30.  and  the  thermal 
versus  mechanical  analogy  for  thermal  conduction,  we 
can  write 

7,  =  }.ca(n.no^.  (33) 

Parameter  p.  which  enters  into  the  function  a  (fl.flo) 
( from  eq  3 1 ).  depends  on  the  hardnesses  of  the  matenal 


sgeggs  C&gy&v  eg  asL  19S9>.!inB^fn»nP=lEo3 
=l3fccso^l^3e^:^scai^ssshss&sm.tof^~4^ofb^ 
tagE^Eqiiafengtg«iiflacoofcniKactegsggil2^^ 
teteiotgi^  besw££a  FCHJ]^)  sad  As 

foiicrasifnaaElaeaSxn'cdiSCTiiwdEbggMl-e^^ 
aEsttegy.  is&  poeo^y  IL  E&e  bc!&  poso^'  He.  n-f^ 
dtsacsesRSS  th?  |»mxn|xsc!ioa  sisortiEie.  sad  ths  pa- 
Eaaa^PafecriaciafcgsfcsaEigyOCEiasfftaifrft  19S6) 
bet«eea^nKdt3:dc^z^iiaiiitzlp{opauesofpo- 
axs.coinsee^ed.pOK^fedtxxais.Tt^coesxiiha- 
nsai  coaducihiy  7^  is  &iammed  esafKikally.  This 
ssdog)*  sid  e!ae  fosnuda  in  eq  33  aie  Ikrdied  by  the 
co£!di3oailmriae3reaoffaggf-{KtiMckccwBar£5besinall 
incoQq]ansoaT«iii3Ebesiz£sof^te£ks.Hence.pon>5ily 
^KNdd  be  in  eIk  nsrge  of  n  ^  0:25  ±  OjQS. 


EXPERDIENTALINVESTIGATIOXS 

Xbo-sicady-siaie  meibods  aie  generally  used  for 
roeasuriae  effeab'c  rhennal  conduaivity  and  thenna! 
di6usmt>’ of  dispersed  and  p(KDUs  media,  as  in  Ar^k- 
pae\’  (196S)  and  Lyalikov  (1965),  for  enample.  These 
meihods  can  be  conveniionally  divided  into  ses'cral 
groups:  1)  meihods  using  mcmoionic  (Plauinov  1973) 
andpenodical  (plippov  19S4)  heating.  2)  methods  us¬ 
ing  equiiiteation  at  an  expcHiential  rate,  e.g..  the  “tegu¬ 
lar  regime  of  Kondratiev”  (Chudnovskit  1962)  and  3) 
probe  meihods  (see.  e.g.,SigaIova  1963.  Bakenov  etaL 
1972,  Zaiuev  ci  al.  1989).  Steady-state  meihods  are 
more  laboiious  and  inaccurate  and  are  therefore  used 
less  often.  For  temperature  intervals  in  which  HM  or 
their  separate  components  (forexample,  thawing  soils) 
undergo  phase  transformations,  traditional  methods  do 
not  work  because  the  thermal  conductivity  equation  is 
nonlinear.  Forsuchcasesihemethodsbasedoncomput- 
er  solutions  of  the  conresponding  inverse  problems  are 
suggested  in  Pavlov  cl  al.  (1980),  Kovalenko  (1986) 
and  Gruzdev  and  Kovalenko  (1988b). 

Experimental  investigations  have,  as  a  rule,  been 
individual,  ad-hoc  efforts.  Their  purpose  in  each  case 
has  been  to  study  thermal  conductivity  of  a  narrow  class 
of  materials.  There  are  practically  no  experimental  in¬ 
vestigations  for  modeled  media. 

Thethcrmal  conductivity  of  soils  (clays,  sandstones) 
and  various  rock  deposits  is  investigated  in  the  works  by 
Bogomolov  (1941),  Sigalova  (1965),  Bakenov  et  al. 
(1972),  Sidorov  (1979),  Nikiolaev  et  al.  (1987)  and 
Zaitzev  et  al.  ( 1 989).  An  attempt  to  summarize  thermal 
conductivity  of  three-  or  four-component  soils  (hard 
frame  with  air,  water  and  oil)  is  made  in  the  works  of 
Volkov  et  al.  (1982).  The  works  of  Franchuk  (1941), 
Kaufman  (1955),  Vasiliev  and  Fraiman  (1967), 
2^brodskii  et  al.  (1968),  Gamashevich  (1974),  Frant- 


sevirii  (1976).  Liansku  mtd  PtsdAeSmch  (I9S2)  and 
Sadov  ( !  9S2)  are  dedicaied  to  experimenta!  study'  of 
the  thennal  condociiviiks  of  theimal  ir»uh!kxi  and 
oebalK3idis^ta3^mb.2askisipk(l9^0)sw!sxaiz.- 
es  experimental  data  for  the  thermal  conductivity  of 
such  systems.  Metal-^etamk  and  ooit^rressed  metal- 
po&tkrmserialsareconsideredinilgewoiksofLyalik- 
<yv  (1965).  SkKokhod  (1967),  Demiddienko  (1972). 
Andre£V’(1975),Aleksaadiovcta!.(i9S5)andGn!Zdcv 
etaL  (1989). 

Experimental  im'cstigatiosis  can  be  unsaiisfacton* 
when  the  materials  studied  are  insufficiently  charaacr- 
ized.  Also,  the  amount  of  experimental  data  is  insuffi- 
deni  to  fav'cr  one  or  another  calculational  model  for 
speciBc  classes  of  porous,  drqreised  materials. 

COXCnLUSION 

Aneview'ofSovietsludiesof thermal  conduaivity  of 
porous  materials  and  soilsshows  that  theoretical  inves- 
tigaiimis  are  more  common  than  experimental.  Many 
models  are  not  adequately  validated  with  the  r.  eager 
experimental  data  that  exist  for  different  classes  of  sys¬ 
tems.  In  faa.  thete  are  no  systematic  experimental  in¬ 
vestigations  of  the  whole  classes  of  HM,  as  they  are 
modeled.  Investigators  give  preference  to  calculations. 
The  lack  of  1)  thermal  measurement  apparatuses,  2) 
accurate  methods  to  characterize  structure,  and  3)  a 
common  data  bank  for  individual  experimental  results 
reinforces  this  syndrome. 

Methods  available  for  determining  the  effective 
ihemtal  conductivity  of  HM  by  theoretical  calculation 
probably  exceed  in  scope  the  whole  range  of  various 
HM  classes.  Choosing  between  models  for  specific 
cases  is  therefore  a  major  problem  in  predicting  of  HM 
ihennal  propeities. 
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